Effects of Noise in Symmetric Two-Species Competition 
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We have analyzed the interplay between noise and periodic modulations in a classical Lotka-Volterra 
model of two-species competition. We have found that the consideration of noise changes drastically 
the behavior of the system and leads to new situations which have no counterpart in the determin- 
istic case. Among others, noise is responsible for temporal oscillations, spatial patterns, and the 
enhancement of the response of the system via stochastic resonance. 
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In the last years the study of deterministic mathemat- 
ical models of ecosystems has clearly revealed a large va- 
riety of phenomena, ranging from deterministic chaos to 
the presence of a spatial organization Jl]-§|]. In the ab- 
sence of spatial dependence these models study the time 
evolution of the averaged number of individuals of some 
interacting species. When space is considered, the de- 
scription is usually done in terms of population density 
fields. These models, however, do not account for the 
effects of noise despite it is always present in actual pop- 
ulation dynamics and arises from different sources, such 
as the intrinsic stochasticity associated to the dynamics 
of the individuals and the random variability of the envi- 
ronment || . Frequently, its effects have been assumed to 
be only a source of disorder ||. Conversely, in this Let- 
ter we show how the consideration of noise in a classical 
Lotka-Voltcrra model of two species competition changes 
drastically and in an unexpected fashion the dynamics of 
the deterministic case. Moreover, consideration of noise 
and space together leads to the appearance of spatiotcm- 
poral patterns which in the deterministic model, except 
for an initial transient and no matter the value of the 
parameters, always looks homogeneous. 

Let us consider a classical Lotka-Volterra model of 
symmetric two-species competition J7j,|| with the addi- 
tion of noise terms 0] defined by the equations 

dx 

— = [j&(l-x- /3y) + f x (x,y)£ x (t) , (1) 

^=Hy(l-y-0x)+fy(x,y)t y (t) , (2) 

where x and y are the population densities, \x is propor- 
tional to the growth rate, and /3 accounts for the inter- 
actions among the species. Here the terms fi{x,y)£i(t) 
(i = x,y) model the contribution of the random forces. 
For the sake of simplicity we assume that f x (x,y) = x, 
f y (x,y) = y. Moreover, is assumed Gaussian 

white noise with zero mean and correlation function 
(€i{i)€j(t + T)) = crS(T)S hj = x.y). This explicit 
form of the noise term may represent, for instance, a 
fluctuating growth rate. Without the noise terms this 
Lotka-Volterra model has been widely studied and it is 
well know that, for (3 < 1, both species are present but 



for > 1, exclusion takes place through a symmetry- 
breaking bifurcation and one of the species is eliminated. 

As an additional feature, we account for a non- 
stationary environment. The simplest way to do it is to 
assume that our dynamics is modulated through a peri- 
odic variation in one of the relevant parameters. Season- 
ality and daily variations are the most important sources 
of regular forcing. The explicit case we consider here is 
described by 

f3 = (3{t) = 1 + e + a sin(w t) , (3) 

where s, a and too are constants. This form of (3(t) corre- 
sponds to a situation in which the competition between 
species is altered in such a way that we move from coex- 
istence (f3 < 1) to exclusion ((3 > 1) in a regular fashion. 
The change in the competition rate may occur when a 
limiting resource which is used by both populations goes 
from large to low values. Thus, at high levels, competi- 
tion is weak, but can become strong at low levels. 

Although some Lotka-Volterra models with sinusoidal 
perturbations of parameters show very complex dynami- 
cal patterns ]8j , in the model we consider the situation is 
rather simple. In the absence of noise, except for a pos- 
sible initial transient, when e > exclusion takes place. 
On the other hand, when e < the species always co- 
exist, in spite of the fact that j3 becomes larger than 1 
periodically in time. The introduction of noise drastically 
changes the previous scenario since even small amounts 
of stochasticity are able to destroy the state of coexis- 
tence periodically |9) . In order to analyze our system we 
have numerically integrated the previous equations by 
using standard methods for stochastic differential equa- 
tions |1(J. In Fig. [j]a we have plotted the temporal evo- 
lution of both species when e < and the noise level 
a is very small, but its effects are slightly appreciable. 
Thus, the species practically coexist all the time. When 
noise increases sufficiently a species is able to nearly ex- 
clude the other one periodically (see Fig. |l|b), i.e., one 
species dominates the other. It is interesting to realize 
that by increasing even more the noise level this coher- 
ent response to environmental variations is lost (See Figs, 
lc and Id). A nontrivial consequence of these results is 
that time variations in competing populations could be 
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FIG. 1. (a) Time evolution of both population densities 
[Eqs. (jjj) and (g)]. The values of the parameters are fj, = 1, 
iuo/2n = 1(T 3 , a = 0.05, s = -0.01, and cr = 10 -12 . (b) 
Same situation as in case (a) but a = 10 -6 . (c) and (d) Same 
situation as in case (a) but a = 10 -2 . 



a noise- induced phenomenon. In this context, periodic 
population oscillations, not present in the original deter- 
ministic model, are now allowed to appear. 

The previous figure exhibits an aspect which has 
received a considerable attention in the recent years; 
namely, the response of the system to a periodic force 
may be enhanced by the presence of noise. This is the 
main characteristic of stochastic resonance (SR) [jy] [l3|| . 
In essence, SR is a nonlinear cooperative effect in which a 
weak periodic stimulus entrains large-scale environmen- 
tal changes. These changes are done in a coherent fashion 
and, as a result, the periodic component is greatly en- 
hanced. The usual form to make the presence of SR man- 
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ifest is through the signal-to-noise ratio (SNR) |11 
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FIG. 2. SNR corresponding to Eqs. (|l|) and (g). The 
values of the parameters are fi = 1, loq/2-k = 10~ 3 , a = 0.05, 
and e = -0.01 






a given quantity describing the state of the system. Thus, 



FIG. 3. Representation of x"j given through Eqs. 
(|) for n = 1750 (a), n = 1875 (b), n = 2000 (c), n = 2125 
(d), n = 2250 (e), n = 2375 (f), n = 2500 (g), and n = 2625 
(h). The values of the parameters are fi = 2, D = 0.05, 
uj /2h = 10~ 3 , a = 0.1, £ = -0.01, and a = 10" 8 . The 
system size is 200 x 200. Black and white colors stand for 
minimum and maximum values, respectively. 



SR arises when the SNR has a maximum as a function 
of the noise level. 

In the system we consider the most straightforward 
form of measuring qualitative changes is by using a quan- 
tity that accounts for the degree of coexistence, such as 
the squared difference of population densities [x — y) 2 . 
In Fig. |^ we have shown the SNR for this quantity which 
clearly exhibits a maximum at non-zero noise level, thus, 
indicating the presence of SR |TJ]. In fact, this is the 
first example of SR in a population dynamics model. In 
biology these situations have been restricted primarily to 
physiological systems but other areas, as popula- 
tion dynamics, have not been explored up to now. We 
should note that previously studied systems displaying 
SR always exhibit oscillations when noise is absent by 
only changing the values of the parameters, usually by 
increasing the amplitude of the input signal. Conversely, 
in the case we are concerned oscillations only appears 
when noise is present. 
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FIG. 4. (a) Representation of x"j given through Eqs. (^) 
and (|) for n = 250, 1250, 2250, 3250 (from left to right). The 
values of the parameters are \i = 2, D — 0.05, ujo/2tv — 10~ 3 , 



a = 0.1, e = —0.01, and a — 10~ . The system size is 
200 x 200. Black and white colors stand for minimum and 
maximum values, respectively, (b) Same situation as in case 
(a) but a — 10 -8 . (c) Same situation as in case (a) but a — 
and e = 0.01. 

A further step in our study will be the analysis of the 
effects of space. The usual way to do it is by adding diffu- 
sive terms to Eqs. (|l|) and (Q). This class of spatiotem- 
poral model corresponds to the situation in which the 
dynamics of the species is continuous in time. Other sit- 
uations of interest concerns with the description of pop- 
ulations whose generations do not overlap in time p7| . 
In this case, the continuous time-space description is no 
longer valid and discrete time evolution models must be 
considered. The usual way to model spatially distributed 
systems whose time evolution is discrete is by using a 
coupled map lattice (CML) |l^|. Here we follow this ap- 
proach. In such a situation the dynamics of our discrete 
model is 

^+ 1 = K,-( 1 -^-/ 3 Xi) 

7 

+ D J2W-y?,j) , (5) 

7 

where (3 n = 1+e+a cos(u;on), D is a constant accounting 
for the diffusion, and indicates sum over the four 
nearest neighbors. Here the random terms are modeled 
by independent Gaussian random variables, denoted by 
XP, and K™, with zero mean and variance unit. The 
remaining parameters have the same meaning as in the 
zero dimensional time-continuous model [Eqs. (Q) and 
(Q)]. When (3 n = this CML model reduces to a logistic 



lattice, which has been widely explored ]19| . 

As in the previous case, when f3 n < 1 this model ex- 
hibits a state in which both species coexist. The co- 
existence state is responsible for the appearance of an 
homogeneous spatial distribution. However, when a suf- 
ficiently large amount of noise is present and (3 n > 1 spa- 
tial patterns arise. In order to elucidate how this spatial 
structure emerges, we have depicted in Fig. || the spa- 
tiotemporal evolution corresponding to the CML model 
for one period of /3". The other periods also exhibit the 
same characteristics. Noise is responsible for the peri- 
odically appearance of the spatial structure since if low 
enough, except for an initial transient, the system always 
looks homogeneous. An example of how patterns are in- 
fluenced by the noise intensity is shown in Figs. fl|a and 
||b. We have plotted four spatial patterns correspond- 
ing to the first four periods of j3 n for two representative 
values of the noise level. The first pattern is strongly 
influenced by the initial conditions, which are random, 
and looks very similar for each noise intensity. However, 
once the initial transient is lost, the patterns correspond- 
ing to the higher noise level are more pronounced than 
the ones with lower noise level. If noise is sufficiently 
decreased patterns do not appear. A similar behavior is 
also present in the continuous time-space model. 

It is worth emphasizing that the deterministic counter- 
part of most models exhibiting noise induced structures 
is able to display similar patterns to those induced by 
noise for a certain range of the values of the parame- 
ters, as for instance the Swift-Hohenberg equation f p0[ . 
In our model, however, these patterns only arise when 
noise is present. For instance if e > a spatial structure 
emerges periodically over a long transient, but eventually 
it disappears. In Fig. we have depicted four patterns 
corresponding to this situation. This figure clearly shows 
how the domains grow in each period. At suffiently large 
time there is one domain, thereby one species excludes 
the other. 

In summary, we have shown that noise cannot system- 
atically be neglected in models of population dynamics. 
Its presence is responsible for the generation of temporal 
oscillations and for the appearance of spatial patterns. In 
contrast, these features do not arise when noise is absent. 
Under some circumstances noise has a constructive role, 
since it is responsible for the enhancement of the response 
of the system via stochastic resonance. The similarity of 
the models we have considered with phase separation, has 
been already pointed out in the literature pl|-p3"| . In this 
regard, it is worth emphasizing that there exist numerous 
systems which can be described through competitive or 
cooperative interactions. To mention just a few: biologi- 
cal assemblies of individuals, coupled chemical reactions, 
political parties, business, and countries |5). Thus, our 
results are not restricted only to population dynamics, 
but the main ideas can be applied to a wide variety of 
situations embracing different scientific areas. 
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